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ABSTRACT 

In this study we compile for the first time comprehensive data sets of solar and 
stellar flare parameters, including flare peak temperatures Tp, flare peak volume 
emission measures EMp, and flare durations r/ from both solar and stellar data, 
as well as flare length scales L from solar data. Key results are that both the 
solar and stellar data are consistent with a common scaling law of EMp oc T^'^, 
but the stellar flares exhibit ~ 250 times higher emission measures (at the same 
flare peak temperature). For solar flares we observe also systematic trends for 
the flare length scale L{Tp) oc T^^ and the flare duration TpiTp) oc T^-^ as a 
function of the flare peak temperature. Using the theoretical RTV scaling law 
and the fractal volume scaling observed for solar flares, i.e., V{L) oc L^'^, we 
predict a scaling law of EMp oc T^^, which is consistent with observations, 
and a scaling law for electron densities in flare loops, rip oc T^/L oc T^-^ . The 
predicted ranges of electron densities are Up ^ cm~^ for solar nanoflares 

at Tp = 1 MK, Up 10^°-" cm-3 for typical solar flares at Tp = 10 MK, and 
Up 10"-^2 Ysjge stellar flares at Tp = 100 MK. The RTV-predicted 

electron densities were also found to be consistent with densities inferred from 
total emission measures, Up = EMp/ qyV , using volume filling factors of qy = 
0.03 — 0.08 constrained by fractal dimensions measured in solar flares. Solar 
and stellar flares are expected to have similar electron densities for equal flare 
peak temperatures Tp, but the higher emission measures of detected stellar flares 
most likely represents a selection bias of larger flare volumes and higher volume 
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filling factors, due to low detector sensitivity at higher temperatures. Our results 
affect also the determination of radiative and conductive cooling times, thermal 
energies, and frequency distributions of solar and stellar flare energies. 

Subject headings: Stellar Flares — Soft X-rays - EUV 

1. INTRODUCTION 

Scaling laws provide important diagnostics and predictions for specific physical models 
of nonlinear processes such as self-organized criticality, turbulence, diffusion, plasma heating, 
and particle accleration. These models have been widely applied in plasma physics, astro- 
physics, geophysics, and the biological sciences. Here we investigate scaling laws of physical 
parameters in solar and stellar flares, which should allow us to decide whether solar and 
stellar flare data are consistent with the same physical flare process. 

The scaling of solar and stellar flare data has been pioneered by Stern (1992), Feldman 
et al. (1995b), and Shibata & Yokoyama (1999; 2002), who showed evidence for a nonlinear 
scaling between the flare volume emission measure EMp and the flare peak temperature Tp. 
These parameters have been measured in solar flares with instruments like Skylab, GOES, 
Yohkoh/SXT (Soft X-ray Telescope) , and RHESSI (Ramaty High Energy Solar Spectroscopic 
Imager), and in stellar flares with ASCA, BeppoSAX, Einstein, EUVE, EXOSAT, Ginga, 
HE AO, ROSAT, Chandra, and XMM-Newton. Compilations of solar flare parameters have 
been presented in Aschwanden (1999), while stellar flare parameters were compiled in a recent 
review by Giidel (2004). In this paper we present for the first time this host of mostly new 
measurements "on the same page" and investigate commonalities and differences between 
the scaling of solar and stellar flares. 

In Section 2 we present the statistical correlations found in stellar flare data, while 
the corresponding counterparts of solar flare data are shown in Section 3. In Section 4 we 
present theoretical modeling of the data, using the well-known RTV law, the generalization 
with gravitational stratification and spatially non-uniform heating, the fractal flare volume 
scahng, and volume filling factor. In Section 5 we discuss the differences between solar 
and stellar scaling laws, the consistency between two different electron density measurement 
methods, and a previously derived "universal scahng law" for solar and stellar flares. Section 
6 summarizes our conclusions. 
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2. STELLAR FLARE OBSERVATIONS 

2.1. Definitions of Physical Parameters 

A stellar coronal flare is usually detected from light curves in extreme ultra-violet (EUV) 
or soft X-ray wavelengths, from which a (background-subtracted) peak count rate Cp [cts s~^] 
at the flare peak time t — tp can be measured. The count rate c{t) for optically-thin emission 

(as it is the case in EUV and soft X-rays) is generally defined by the temperature integral of 
the total (volume- integrated) differential emission measure distribution dEM{T)/dT [cm~^] 
and the instrumental response function R{T) (in units of [cts s~^ cm^]), 

i-d- c{t) = / '^^^ RiT) dT . (1) 

where the factor {And'^) comes from the total emission over the full celestial sphere at a 
stellar distance d (in parsecs). The differential emission measure distribution (DEM) of 
flares shows usually a single peak at the flare peak temperature Tp, so that the emission 
measure peak at the flare peak time, EMp = dEM{t = tp,T)/dT ^ dEM{t = tp,T = Tp), 
can be approximated with a single temperature (which corresponds to an emission measure- 
weighted average value), 

And'^ Cp = 47rd^ c{t = tp) ^ EMp R{Tp) . (2) 

The total (volume-integrated) emission measure EMp at the flare peak is deflned as the 
squared electron density n integrated over the source volume V, 

EMp = Jn^dV^nlV , (3) 

where the right-hand approximation implies that = n^(t = tp,T = Tp) is the squared 
electron density at the flare peak time averaged over the volume V of the flare plasma, 
assuming a unity fllling factor. 

Integrating the count rate c{t) over the flare duration rj yields the total counts C, which 
in the case of a single-peaked DEM can also be approximated (with Eq. 2) as 

4rfC = 4rf/cW*.4rfc,.,^^M,i^ra.,. (4) 

The radiative loss rate for optically thin plasmas is a function of the squared density 
and the radiative loss function A(T), 

^ = n^nMT) « , (5) 
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(in the coronal approximation of fully ionized plasma, i.e., Ug ~ rii) where the radiative loss 
function has a typically value of A(T') io~23...-22 ^^^^ g-ij ^^le temperature range 
of T f» 10^ ••8 K. From this we can define a peak luminosity Lx in soft X-rays by integrating 
over the volume and temperature range, 

Lx = V J n'^{t = tp, T)K{T)dT « EMp A{Tp) . (6) 

The total radiated energy Ex integrated over the flare duration is then 

Ex^ J J J n^{t,T)A{T) dV dT dt « EMp A{Tp) Tf . (7) 

This yields a convenient conversion from observed total counts iird'^ C (Eq. 4) into total 
radiated energy Ex (Eq. 7), 

K[lp) 

which involves a temperature-dependent conversion factor f{Tp) — A{Tp)/R{Tp). 

For comparison we calculate also the total thermal energy Et of the flare volume at the 
flare peak time t — tp, 

/Q / TT' 1\ /T T' 
3n{t = tp, T)kBT{t = tp)V{t = tp)dT ^ SnpheTpV = — ^ (9) 
Up 

where Up = n(t = tp,T = Tp) represents the electron density at the flare peak time t = tp 
and DEM peak temperature T = Tp. The relation between the total thermal energy E^ and 
the total radiated energy Ex is then 

Et - Ex——%^^--— , (10) 
np{Tp) A{Tp) Tf[Tp) 

where the peak electron density np{T^ and the flare duration Tf{Tp) may have a statis- 
tical dependence on the flare peak temperature Tp, and this way define the temperature 
dependence in the correlation between the thermal energy Et and the total radiated energy 
Ex- 



2.2. Statistical Correlations 



A recent compilation of stellar flare measurements is given in Giidel (2004; Table 4 
therein). This database contains measurements of the total energy radiated in soft X-rays, 
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the peak emission measure EMp, flare peak temperature Tp, and flare duration t/ 
from 68 different stellar flares. Generally, only the peak emission measure EMp (derived 
from the peak count rate Cp [Eq. 2] and the known response function R(Tp)), the flare peak 
temperature Tp (if a EUV/soft X-ray spectrum is available), and the flare duration Tf can 
be measured directly, which we will consider as independent variables in the following, while 
all other quantities such as the energies {Ex,Et) and densities Up are derived quantities, 
using relations as given in §2.1. 

We show scatterplots of various parameters as a function of the flare peak temperature 
in Fig. 1. Wc perform linear regression flts of y{x), x{y) (shown with thin lines in Fig. 1), with 
the ordinary least square bisector method (shown with thick lines in Fig. 1), and calculate 
the linear regression coefficients (indicated with RC in Fig. 1). We find a highly significant 
correlation (with a regression coefficient of RC — 0.68) for the fiare peak emission measure 
EMp as a function of the fiare peak temperature Tp (Fig. 1 top left), i.e., 

EM,(T,) = 10^°-« (^^^) , « = 4.5 ±0.4. (11) 

This value is essentially identical to the correlation given in Giidel (2004), where a powerlaw 
slope of q; = 4.30 ± 0.35 is quoted, determined with a similar hnear regression method. We 
test this general EMp — Tp scaling law in evolutionary curves observed from 8 different stars 
in the next Section (§2.3). 

The other independently measured parameter is the fiare duration t/, for which we find 
a marginally significant correlation {RC — 0.39) with the fiare temperature Tp (Fig.l lower 
left), 

^^^^-) = '°''(iOMk)' ' = (12) 

The other correlations can be understood as a consequence of the correlations found 
between the independent parameters Tp, EMp, and Tf. For instance for the peak X-ray 
luminosity Lx we find (Fig. 1 top right), 

Lx(r,) = 10^^« , 7-4.7±0.4. (13) 

because the peak luminosity (Eq. 6) is proportional to the peak emission measure, i.e., 
Lx oc EMp, if we neglect the weak temperature dependence of the radiative loss function 
A(rp), and thus should have about the same powerlaw slope, i.e., 7 = 4.7 ^ a — 4.5. 

For the total radiated X-ray energy Ex we find a correlation of (Fig. 1 bottom right), 
^x(r,) = 10^°-^ , 5 = 6.1 ±0.5. (14) 
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which is expected to scale as Ex{Tp) oc EMp{Tp)Tf{Tp) oc (Eq. 7) with 5 ^ a + (3 ^ 
4.5 + 1.8 = 6.3, which indeed agrees with the best fit, 6 — 6.1 ± 0.5. 

We cannot calculate the thermal fiare energy Et (Eq. 10) with current data, since 
we need the knowledge of the fiare peak density Up, which is generally not independently 
measured from EMp and Tp. The values quoted in Table 4 of Giidel (2004) have a mean 
value of Up = io^^-^='=*'-^ cm^'^ and show no temperature dependence. The fiare peak density 
Up can be calculated with relation (Eq. 10), if one uses some theoretical assumptions (e.g., 
RTV scaling law, see §4.1) which predicts the temperature dependence of the electron density 
np{Tp). 



2.3. The Evolution of Stellar Flares 

In the EMp — Tp scaling law (Eq. 11), the observables represent the peak emission 
measure EMp = EM{t = tp,T = Tp) and peak temperature Tp = T(t = tp) measured at the 
peak time t = tp of stellar flares. However, since both observables EM{t) and T{t) change 
during the flare as a function of time, we inquire how closely the temperature predicted 
by the statistical scaling law (Eq. 11), 

(EM \ '^'^^ 
-^J , To^lOMK, EMo^W' cm-\ (15) 

matches the observed peak temperature Tp at the peak time t = tp, and how the peak 
temperature Tp and the statistically predicted peak temperature T5 matches the maximum 
temperature Tm of observed stellar flares. For this purpose we collected evolutionary phase 
diagrams of the temperature T{t) versus the emission measure EM{t) observed in 8 different 
flares (Fig. 2). The curves of 5 stellar fiares observed with GINGA are from the stars HR1099 
[or V711 Tau] (Stern 1996), 1EQ1839.6+8002 (Pan et al. 1997), II Peg (Doyle et al. 1991), 
UX Ari (Tsuru et al. 1989), and Algol (Stern et al. 1992); two observations of AB Dor 
fiares are observed with BeppoSAX (Maggio et al. 2000), and one observation of Proxima 
Centauri was obtained with XMM-Newton (Reale et al. 2003). The evolutionary curves 
are shown in Fig. 2, including error bars in temperature, and some are shown with 90% 
confidence regions for both derived (EM and T) parameters. The peak emission measures 
EMp are indicated with a vertical dashed line in Fig. 2, ranging from EMp = 10^^'^ cm~^ 
(for Proxima Centauri) up to EMp = 10^^'^ cm^^ (for a large fiare on UX Ari). The rise 
time has not been captured for two observations (II Peg and UX Ari). The observed peak 
temperatures Tp are also indicated, ranging from Tp — 10^-^° = 12.5 MK (for Proxima 
Centauri) to Tp = 10^ °^ = 115 MK (AB Dor, 1997-Nov-29). At the same time we plot also 
the statistically predicted scaling law temperatures Ts (using Eq. 15) and the fiare maximum 
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temperatures Tm, which is usually reached shortly before the emission measure peak. Note 
that the flare rise covers the upper (high-temperature) part of the evolutionary curve, the 
flare peak time t = tp is defined at the right-most datapoint (highest emission measure), 
and the flare decay or cooling phase extends over the lower part of the evolutionary curve, 
from right to left, so the time t can be tracked in clock-wise direction along the evolutionary 
curve. We list the three values of temperature Tp, T^, and Tm in Table 1 and see that there 
is a close match between them, differing not more than \log{Ts/Tp)\ < 0.13, corresponding 
to a factor of ^ 1.3. This means that the effective flare peak temperature Tp is predicted 
with an uncertainty of ~ 30% using the observed flare peak emission measure EMp and the 
statistical EMp — Tp scaling law (Eq. 11). 



3. SOLAR FLARE OBSERVATIONS 

3.1. Flcire Peeik; Emission Measure versus Temperature 

In Fig. 3 we show a compilation of flare peak emission measures EMp versus flare peak 
temperatures Tp observed in solar flares. We included data sets from small-scale flares (also 
called nanoflares) in EUV to large-scale flares observed in soft X-rays. Small-scale flares (also 
called heating events in quiet corona) were measured in EUV with typical temperatures of 
T ^ 1.0 — 1.5 MK (derived from 171/195 A fllter ratios) and total emission measures of 
EMp ^ 10"^^ - 10"^^ cm-3 (Krucker & Benz 2000; Aschwanden et al. 2000a). The next 
larger category of small flares is observed in soft X-rays, also called active region transient 
brightenings, for which Shimizu (1995) measured typical peak emission measures of EMp 
1044.5 _ 1047-5 pgg^j^ tcmpcraturcs of Tp 4 - 8 MK with Yohkoh/SXT. For large 

solar flares, peak emission measures in the range of EMp ~ 10^^ — 10™ cm~^ and flare 
peak temperatures in the range of Tp f« 6 — 30 MK were measured, using observations 
from Skylab (Pallavicini, Serio, & Vaiana 1977), GOES (Feldman et al. 1995a, 1996; Garcia 
1998), Yohkoh/BCS (Sterling et al. 1997), and RHESSI (Battagha, Grigis, & Benz 2006). 
Comparison of flare temperatures simultaneously measured with GOES and RHESSI reveal a 
systematic bias that RHESSI determines higher temperatures, i.e., Trhessi ~ 1-31 Tgoes-\- 
3.12 [MK], which indicates that RHESSI Hts only the high-temperature tail (at > 3 keV), 
while the GOES temperatures are weighted by the peak of the emission measure distribution 
(Battaglia et al. 2006). Combining all these measurements together (except for RHESSI, 
which has a particular high-temperature bias) and performing a linear regression flt in the 
entire temperature range of 7], ~ 1 — 30 MK with the ordinary least square bisector method 
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we find a statistical correlation (Fig. 3) of 

EMM = lO'^-' (lo^) " [^"^"'1 , a = 4.7 ± 0.1 . (16) 

Note that the powerlaw slope of the EMp—Tp relation is essentially identical for both the solar 
and stellar fiare data sets, but stellar flares have a higher temperature range [Tp ^ 10 150 
MK) and have about a factor of io^°-*~^^-^ ?a 250 (comparing the factors in Eqs. 11 and 16) 
higher emission measures in the overlapping temperature range {Tp fa 10 — 30 MK). 

3.2. Flare Duration versus Temperature 

In Fig. 4 we show a compilation of flare durations r/ versus flare peak temperatures 
Tp observed in solar flares. The flare durations were obtained from Tf = {tend — tstart) of 
12 TRACE nanoflares (Fig. 9 of Aschwanden et al. 2000a), from Tf = Ttotai / Speaks of 23 
SoHO/EIT quiet-Sun brightening events (Table 1 of Krucker & Benz 2000), from 200 
active region transient brightening events observed with Yohkoh/SXT (Fig. 3 in Shimizu 
1995), from the total durations tt of 31 flares observed with Skylab (Table 3 of Pallavicini 
et al. 1977), from t/ = (r^se + 'Tdecay) of 14 flares observed with Yohkoh/SXT (Table 6 of 
Garcia 1998), from Tf = {tend — tstart) of 9 GOES light curves (Sterling ct al. 1997). and 
from 19 flares observed with Yohkoh/SXT (Table 1 of Mctcalf & Fisher 1996). The EUV 
nanoflares have a typical duration of rj ^ 5 x 10^ — 2 x 10'^ s (~ 7 — 33 min), while the larger 
flares have typical durations of Tf ^ 5 x 10^ — 3 x 10^ s (^10 min - 10 hrs). Solar and stellar 
flares have comparable durations in the overlapping temperature range of T ^ 10 — 30 MK. 
Fitting all solar and stellar flare durations combined with the ordinary least square bisector 
method we obtain a statistical correlation (Fig. 4) of 

rf{Tp) = 10=^-^ (l^) ^ i'h 0.91 ± 0.05 . (17) 

which encompasses 67% of the flare durations within a factor of 3. 

3.3. Flcire Length Scales versus Temperature 

We can measure geometric parameters of flares only in solar data where we have spa- 
tial resolution using EUV and soft X-ray imagers. The most directly measured geometric 
parameter is the length L or area A of a flare, while the flare volume V can only be inferred 
indirectly from the flare area. We compile measurements of the flare length scales L versus 
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flare peak temperature Tp in Fig. 5, where we quote either the observerd (projected) loop 
lengths L, or length scales converted from the measured area A, i.e., L — A^/^, or quoted 
volume, i.e., L oc V^^'^. 

Flare length scales were measured in the range of: L — 2 — 11 Mm for 12 nanoflares 
observed with TRACE {L = ^/lxw in Table 1 of Aschwanden et al. 2000a), L ^ A - U 
Mm for 23 quict-Sun brightenings observed with SoHO/EIT (from L = ^/A in Table 1 of 
Kruckcr & Benz 2000), L = 3 — 36 Mm for 23 soft X-ray bright points observed with MSSTA 
{Multispectral Solar Telescope Array) (half lengths L in Table 1 of Kankelborg et al. 1997), 
L ^ 5—40 Mm in ^ 200 transient soft X-ray brightenings observed with Yohkoh/SXT (Fig. 4 
in Shimizu 1995), L = 5 — 54 Mm in 31 flares observed with Skylab (from L = V^l"^ in Table 
3 of Pallavicini et al. 1977), L = 12 - 144 Mm for 14 flares observed with Yohkoh/SXT 
(observed loop lengths Lq^^ in Table 6 of Garcia 1998), L = 3 — 13 Mm in 20 soft X-ray 
flares observed with Yohkoh/SXT (Table 1 in Reale et al. 1997), and L = 3 - 31 Mm in 19 
soft X-ray flares observed with Yohkoh/SXT (Table 1 in Metcalf & Fisher 1996). Fitting 
all solar length scales combined with the ordinary least square bisector method we obtain a 
statistical correlation (Fig. 5) of 

L{Tp) = 10'-' (y^) ^ [cm] , P = 0.91 ± 0.04 . (18) 
which encompasses 67% of the flare durations within a factor of ~ 3. 

4. THEORETICAL MODELING OF SCALING LAWS 

4.1. The Standard RTV Law 

A scaling law between the peak temperature Tmax, pressure p, and loop length L in 
static coronal loops has been derived by Rosner, Tucker, & Vaiana (1978), assuming an 
equilibrium between a spatially uniform heating rate (Eh) and the conductive (Econd) and 
radiative (Erad) loss rates, the so-called RTV scahng law, 

T^,,{p,L)^im{pLy/\ (19) 

While the validity of this scaling law applies to a coronal loop in hydrostatic equilibrium, it 
might also approximately apply to a flaring loop near the peak time, because both (1) the 
energy and (2) momentum equations are nearly balanced near the flare peak. 

(1) Energy equation: In the initial heating phase, the heating term is larger than the 
combined conductive and radiative loss term, i.e.. Eh > {Econd + E^ad), because the average 
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flare temperature is rising {dT{t) /dt > 0) due to excessive heating. After the peak time, the 
conductive and radiative losses exceed the decreasing or stopped heating rate, i.e.. Eh < 
{Econd + Erad), bccausc the flare plasma is cooling {dT(t)/dt < 0). Therefore, there is a 
balance between heating and loss terms, i.e.. Eh ~ {Econd + Erad), near the flare peak 
time tp [dcflned by the peak in emission measure EMp = EM{t = tp)]. Hydrodynamic 
simulations by Jakimiec et al. (1992) show that the RTV scaling law predicts a maximum 
electron density Up that can be reached in a flare loop if a constant heating rate is applied 
sufficiently long after the maximum temperature Tp is reached. In Appendix A we show 
how the agreement between the maximum electron density Up and the RTV-predicted value 
scales with the heating duration, but is nearly independent of the maximum temperature Tp 
and heating rate Eh- 

(2) Momentum equation: Secondly, also the momentum equation is nearly balanced 
after a flare loop is flUed. The losses are dominated by thermal conduction at high flare 
temperatures (say at T > 10 MK; Aschwandcn & Alexander 2001), and the loop filling time 
is much shorter (in the order of Tfm < 1 min; see hydrodynamic simulations by MacNeice et 
al. 1984; Nagai & Emshe 1984; Fisher et al. 1985a,b,c; Mariska & Poland 1985; Yokoyama & 
Shibata 1998, 2001; Hori et al. 1997, 1998) than the conductive or radiative cooling time (in 
the order of TcooI ~ 10 min; Antiochos & Sturrock 1978; Culhane et al. 1994; Aschwanden & 
Alexander 2001), and thus pressure gradients resulting from the chromospheric evaporation 
process are largely balanced out so that the momentum equation is approximately fulfilled. 
Also the assumption of a constant pressure (made in the RTV law) is better fulfilled in the 
hot soft X-ray emitting flare loops (because of the larger pressure scale heights. At > L) 
than in the cooler EUV-emitting coronal loops (where often Xt < L). 

Applying the standard RTV scahng law to flare loops now, with Tp Ri T^ax, by inserting 
the expression for the thermal pressure, 

p(np, Tp) = 2npkBTp , (20) 

we obtain a scaling law for the peak density Tip QjS Qj function of the peak temperature Tp 
and loop length L, 

np{Tp, L) = Co ^ , Co = Yiooh^ [cm-^R-^] . (21) 

where we defined a constant Cq for the numerical factor. The electron density Up, however, 
is not a directly measured quantity in most solar and stellar flare observations, so instead 
we use the directly measured quantity of the total (volume-integrated) emission meausure 
EMp, 

EMp{np,V)^nlV , (22) 
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involving the total volume V of all flare loops. The standard RTV law, expressed in terms 
of the total emission measure EMp, is then 

EMp{Tp, L, V) = nliTp, L) V ^ 4 ^ . (23) 



4.2. Serio's Scaling Law 

It is known that the RTV law underestimates the densities of active region loops observed 
in EUV (at temperatures of T 1—2 MK), as a result of the (sometimes invalid) assumptions 
of a spatially uniform heating rate and constant pressure (Lenz et al. 1999; Aschwanden, 
Nightingale, & Alexander 2000b; Aschwanden, Schrijvcr, & Alexander 2001; Aschwanden 
et al. 2003; Winebarger et al. 2003; Aschwanden ct al. 2007). Applying Serio's scaling law 
(1981), which generalizes the RTV law by including a non-uniform heating scale height Sh 
and gravitional stratification, leads to the following correction factor qserio, 

QSerio i 

(24) 

qserioiL, Tp, Sh) = exp ( -0.08— - 0.04^ ) . (25) 

where L is the loop half length, sh the heating scale height, and \p the pressure scale height. 
This correction factor depends on (L, Tp, sh), since the pressure scale height is a function of 
the temperature, Ap 47 x Tmk [Mm]. In order to apply Serio's scaling law we can simply 
replace the constant Cq with a new function Ci(L, Tp, sh), 

Ci{L,Tp,SH) = = Co exp ("0.24— + 0.12^) . (26) 

(ISerio \ Sh / 

If we insert Serio's correction factor Ci = Co/glg^jg (Eq. 26) into the RTV law of the total 
emission measure (Eq. 23) we obtain the following scaling, 

EMp{Tp, L, V, sh) = cl = cl'^ exp (0.48— + 0.24^) . (27) 

Li L Sh \ 

Stellar flares have temperatures of Tp > 10 Mm, and thus pressure scale heights of Ap = 47 
Mm x(Tp/l MK) >470 Mm that are much larger than the expected solar loop lengths 
(typically Lq = 25 Mm at To=10 MK). We can therefore assume L <^ Xp and neglect the 
second term exp(L/Ap) in Serio's correction and can approximate (Eq. 27) by 

EMp{Tp, L, V, Sh) ^ c^^Z gxp (0.5—) . (28) 

Sh 
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which yields higher emission measures than the RTV law (Eq. 23) for heating scale heights 
of sh ~ L. However, for relatively short heating scale heights, say Sj^/L < 1/3, a density 
inversion occurs at the temperature maximum of the loop, which is unstable against the 
Rayleigh- Taylor instablity (Vn ■ g < 0). This is confirmed by analytical calculations as 
well as with hydrodynamic simulations, where a lower limit is found for stable hydrostatic 
equilibria at sh/L < 1/3 (Serio et al.l981; Aschwanden & Schrijver 2002; Winebarger et 
al. 2003). Given this instabilitiy limit we find that Serio's correction factor can increase 
the emission measure EMp only by a factor of exp(0.5 x 3) 4, or the electron density 
He oc V EM by a factor of ~ 2, which is also confirmed by hydrodynamic simulations 
(Winebarger et al. 2003). 

4.3. Fractal Scaling of Flare Volume 

In early models, stellar flares were modeled with a single loop, so that the volume 
V — L X a was described by a loop length L and a cross-sectional area a. If this cross- 
sectional area a is chosen as a constant for small and large flares, the flare volume would 
just scale hnearly with the length scale L, i.e., V oc L^. However, we think that such 
geometries are unrealistic, given the fact that large solar flares always reveal arcades with 
hundreds of flare loops, remnants of multiple magnetic reconnection sites. A uniformly fllled 
flare arcade is expected to scale with oc in the Euclidian limit. In arcsccond high- 
resolution images such as from TRACE, however, it becomes evident that the flare region 
is not uniformly filled, but rather has a filamentary structure that can be described with a 
fractal dimension or an area filling factor. A detailed study (Aschwanden & Aschwanden 
2006a) of 20 GOES X- and M-class flares observed with TRACE has revealed that the fractal 
area (normalized by the time- integrated flare area Af) varies from near zero at the beginning 
of the flare to a maximum of A(t)/Af = 0.65 ± 0.12 at the peak time of the flare, which 
corresponds to an area fractal dimension of Da $ 1.89 ± 0.05 (at the flare eak time), also 
called Haussdorf dimension Da, which deflnes the scaling of a fractal area A with length scale 
L, i.e., A{L) oc L^^. Also a statistical study of nanoflares has shown a Haussdorf dimension 
of Da = 1.5 ±0.2 (Aschwanden & Parnell 2002). From the measured area fractal dimensions 
Da at the flare peak time, a volume fractal dimension Dy ~ 2.37 ± 0.14 was derived for a 
flare model with an arcade geometry (Aschwanden & Aschwanden 2006b), where the fractal 
flare volume scahng is deflned by V{L) oc L^^ . The resulting volume flUing factors were 
found to be in the range of qy ~ 0.03 — 0.08 at the flare peak time. 

The implication of the previous studies is that the flare volume has a fractal scahng of 
V{L) oc L^-^, rather than Euclidian flUing with V{L) oc L^, which affects our scahng laws. 
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We can quantify this fractal volume scaling by 

V{L) = qyLl {-] , (29) 



where Lq = lO^ "^ cm (=25,000 km) is the average length scale of a solar flare loop (at a 
typical flare temperature of Tp = 10 MK, according to Eq. 18), qv = 0.03 — 0.08 is the 
average volume filling factor, and Dy ~ 2.37 ± 0.14 is the volume fractal dimension. This 
definition fulfills the normahzation V{L — Lq) — qyL^, so the average volume filling factor 
is defined by qv — V/Vq at Vq — L^. 

Inserting this volume scaling V into the RTV law (Eq. 23) we obtain the following 
relation, 

EMp = cl%^^ clqvUT^ (pi ' oc T^L°-^ . (30) 



Therefore, comparing total emission measures of solar and stellar fiares at the same fiare 
peak temperature, we expect the following dependence 



em; _ q^ /L*^""^-' 



EMf 5® V L® 



(31) 



while the relative densities scale reciprocally to the loop length L (for the same fiare peak 
temperature) according to the RTV law (Eq. 21), 



(32) 



Is this fractal volume model consistent with the scaling laws we found for solar flares? 
Inserting the observed scaling of L{Tp) a; T°-^ (Eq. 18 and Fig. 5) into the relation for the 
total emission measure (Eq. 30) we find the following scaling with the fiare peak temperature 
T 

EMp oc [L(7;)]^--2 oc T^-^ . (33) 

which is close to the observed scahng of EMp{Tp) oc T^-'^^^-'^ for solar fiares and EMp{Tp) oc 
rp4.5±0A fQj. stellar fiares (Fig. 3), and thus the fractal scahng is consistent with observations. 
This RTV model (Eq. 21) predicts also the following dependence of the density Up with fiare 
peak temperature Tp, using the scaling L{Tp) oc T°-^ observed in solar fiares, 

VT,)^c„J|=10-(^)"[cm-l. (34) 

Thus, we predict for stellar fiares in the temperature range of Tp fa 10 — 100 MK up to an 
order of magnitude higher electron densities than for solar fiares in the temperature range 
of Tp^ 10- 20 MK. 
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4.4. Flare Conductive Cooling Times 

For the physical process of conductive coohng, there is a simple relation between the 
maximum temperature Tjnax and loop length L for a coronal loop in hydrostatic equilibrium, 

dtdV^-^'^^- dT^ ds^" 7' ds t'^L^ ' ^^^^ 

where k = 9.2 x 10^^ [erg cm~^ K^^/^] is the classical Spitzer conductivity. We can argue 
that the heating rate is dominated by conductive losses at the flare peak time, if the flare 
is sufficiently hot so that radiative losses can be neglected. For solar flares it has indeed 
been shown that conductive cooling always dominates over radiative losses in the hot phase 
of flares, say at temperatures of T > 10 MK (e.g., Antiochos & Sturrock 1978; Culhane et 
al. 1994; Aschwanden & Alexander 2001). Assuming this balance of heating and conductive 
loss at the flare peak time, i.e., dEu/dt dV ~ dEcond/dt dV, leads (by inserting the first 
RTV law of Eq. 19 into Eq. 35) directly to the second RTV scaling law for the heating rate 
(Rosner et al. 1978), 

dEn dE,ond T'l^ 

d^W -l^tW ^ IN-^ ■ 

The thermal energy of the flare plasma per volume is, 

dE'^ 

— — = SnphsTp , (37) 
dV 

which can be expressed in terms of flare peak temperature Tp and flare length scale L (by 
substituting the flrst RTV law Eq. 21), 

"^^^ Shco (^) . (38) 



dV \L ^ 

The two expressions of the thermal energy (Eq. 38) and the conductive cooling rate (Eq. 35) 
allow us to deflne the cooling time Tc by thermal conduction, 

dE^/dV _ 3hco (T^/L) ^ L 

If we express this relation in dimensionless units in terms of the reference values Tq = 10^ K 
and Lq — 25 Mm we obtain the following scaling law, 

tJT^,L)^t,o (—](—] ~ 1600 ( — ( [s] . (40) 

In Fig. 6 (top left panel) we plot this theoretically estimated cooling time Tc{Tp, L) calculated 
from the observed values of the flare peak temperature Tp and flare length scale L as a 
function of the observed flare duration times t/. We flnd that the theoretically estimated 
conductive cooling time Tc of most flares is comparable or shorter than the flare duration r/. 
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4.5. Flare Radiative Cooling Times 

The radiative loss rate is generally expressed as a product of the electron density ng, ion 
density n^, and the radiative loss function A(T), which in the coronal approximation (with 
full ionization, i.e., rie ~ n,) is 



dV dt 



ne, T) = niK{T) , (41) 



where the radiative loss function can be approximated by piece-wise powerlaws (Rosner et 
al. 1978; Mewe et al. 1985), 

( 10-21-94 for 105-^5 ^rj. ^ -^q6.3 k 

A(T) = I 10-17-73T-2/3 for 106-3 < T < 10^-^ K . (42) 

I 10-24-66j^l/4 foj. J. > 10^-3 K 

The radiative cooling time can then be defined as the ratio of the thermal energy (Eq. 39) 
and the radiative loss rate (Eq. 41), where we can eliminate the unknown density by inserting 
the RTV law (Eq. 21), 

dE^/dV _ 3hco {T^/L) _ 3h L 
"■"^^ ^' ^ dEraa/dtdV clK{T,)(T^/L^) c,K{T^)T^' ^ > 

which reads in dimensionless units, 

' 700 (^) (^) ^ for 10^-^^ <T< 106-3 K 
TriT,, L) = <! 2000 {^) {-^)"'' for 10^-3 < T < 10^ K [s\ . (44) 
6500 (^sii^) (^y^' for T> 10^-3 K 



We plot the radiative coohng times Tr{Tp, L) calculated with the observed values (Tp, L) in 
Fig. 6 (top right panel) as a function of the flare duration t/. These radiative cooling times 
based on the RTV scaling law are generally longer than the conductive cooling times, and 
they clearly exceed the flare durations for most of the flares, up to an order of magnitude for 
the short EUV flares with durations of Tf < 103 s. Since the observed flare duration should 
be an upper limit of the cooling time, either the conductive or radiative cooling time should 
be equal or shorter. A combined cooling time r can be deflned from the exponential folding 
time that would result from the product of the two exponential cooling processes, 

111 

- = - + -. 45 
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We plot this combined cooling time r as a function of the flare duration in Fig. 6 (bottom left 
panel) using the RTV law. The so-defined combined cooling time is almost always shorter 
than the observed flare duration. 

Including Serio's correction factor (Eq. 25) for the conductive and radiative cooling 
times, 

^Serio _ ^RTV ^1 



^Serio ^RTV '^0 
T — T — 

Ci 

yields only small corrections, shown for the five solar datasets in Fig. 6. The average cor- 
rection values are < qserio >= 0.90 — 0.98 for the EUV datasets (Aschwanden et al. 2000a; 
Krucker & Benz 2000), and < qserio >= 0.05 — 0.95 for the soft X-ray datasets (Pallavicini 
et al. 1977; Garcia 1998; Metcalf & Fisher 1996). We plot the corrected flare cooling times 
predicted by Serio's scaling law in Fig. 6 (bottom right panel). The major effect of Serio's 
correction is that the cooling times of the larger and hotter soft X-ray emitting flare loops 
become shorter, limiting essentially all flare loop cooling times to r < 10^ s. Thus, large 
flares that last signiflcantly longer (up to t/ < 2 x 10^ s) must consist of multiple subflares. 



rr""-^ , (46) 



1 

^Serio 

■ = QSerio ■ (47) 



5. DISCUSSION 

5.1. Scaling Laws of Solar versus Stellar Flares 

The combined EM—T diagram for solar and stellar flares (Fig. 3) shows some similarities 
but also intriguing differences in the scaling behavior. Solar ffares have been observed mostly 
in the temperature range of Tp ^ 7 - 30 MK (or Tp ^ 1 - 30 MK if we include the EUV 
nanoflarcs), while stellar flares have been detected within the temperature range of Tp ^ 
10 — 150 MK. The lack of observations of cooler stellar flares Tp < 10 MK is likely to be due 
to the sensitivity limit, which is about at EMp > 10^^ cm^^ for stellar flares. The sensitivity 
limit also systematically increases with higher flare temperatures, up to ~ 10^^ cm~^ for 
the hottest stellar flares with Tp > 100 MK, which is hkely a consequence of the decreasing 
sensitivity of current soft X-ray detectors at higher temperature lines. 

What is similar for both solar and stellar flares is the overall slope of the EMp — Tp 
relationship, which was found to have a powerlaw slope of a = 4.7 ±0.1 for solar flares 
(excluding the RHESSI data that have a high-temperature bias), and a slope of slope of 
a = 4.5 ± 0.4 for stellar flares (Fig. 3). Both powerlaw slopes are consistent with the 
theoretically expected scaling of a 4.3 (Eq. 33), which is based on the RTV scaling law 
(Eq. 23), the fractal volume scaling (Eq. 29), and the observed spatial scaling L{Tp) T°-^ in 
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solar flares (Eq. 18). This agreement supports the assumption of the fractal volume scaling. 

If we were to assume a monolithic single-loop model with constant cross-section that scales 
with V{L) oc L, the resulting scaling law would be EMp oc T^'^, or a monolithic cubic 
model with V{L) oc would yield EMp oc T^^, which are both less consistent with the 
observations. 

What is different between solar and stellar flares is the emission measure in the same 
temperature range. In the overlapping temperature range of Tp ^ 10 — 30 MK we find that 
the total emission measure of stellar flares is larger by an average factor of 

EM* 1050-8 

bJsI = 1551-250, (48) 

comparing the factors (Eqs. 15, 16) of the two linear regression fits in Fig. 3. If we compare 
the theoretical model of the RTV law and the fractal volume scaling (Eq. 31) for the same 
temperature range (Eq. 33), we see that the total volume- integrated emission measure ratio 
of stellar to solar flare depends on the volume fiUing factor gy, length scale L, and volume 
fractal dimension Dy/. Since the sensitivity limit of stellar soft X-ray detectors represents a 
bias for larger emission measures, we think that the observed stellar flares have a bias for 
both higher volume filling factors qv and flare sizes L. For solar flares we found spatial filling 
factors in the range of qv = 0.03 — 0.08 (Aschwanden & Aschwanden 2006b), so stellar filling 
factors can be up to a factor of qy/qy — 30 higher, ameliorating the size requirement to 
be (L*/L®)'^^~^ > 10 — 25 to match the 250 (Eq. 48) times larger emission measures of of 
detected stellar flares . 



5.2. The Electron Density in Solar and Stellar Flares 

The electron density Hp of the flare plasma cannot directly be measured and is therefore 
dependent on the volume model if derived from the total emission measure. If the density 
is naively derived from an Euclidian flare volume, i.e., Up = y^EMp/V, we obtain only a 
lower limit. However, detailed measurements of the area fractal dimension and modeling of 
the volume fractal dimension has yielded volume filling factors of qy ~ 0.03 — 0.08 for solar 
flares, which raises the average electron densities in the flare loops (at the flare peak time) 
by a factor of 1/^/gy ~ 4 — 6. Thus, a more realistic estimate of the plasma density in flare 
loops, based on the observed total emission measure EMp and flare length scale L is (using 
the fractal scaling of Eq. 29) 

"^■^^-( ..LgfvV )"'' '''' 
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where the volume filling factor is typically in the range of qv ~ 0.03 — 0.08 and the fractal 
dimension is Dy ^ 2.4 according to measurements of fractal flare areas (Aschwanden & 
Aschwanden 2006a,b). 

On the other hand we can estimate the plasma density theoretically, using the RTV law 
(Eq. 21) apphed to the observed flare peak temperature Tp and loop half length L, 

nf''{Tp,L)^Co^. (50) 

We show the two estimated electron densities versus each other in Fig. 7 and find a good 
agreement within less than an order of magnitude. Both methods yield flare densities in the 
range of rie ~ 10^ — 10^^ cm~^ for solar flares. 

Plugging in the observed scaling law for loop lengths versus flare peak temperature, 
i.e., L{Tp) oc Tp-^ (Eq. 18), we find an approximate prediction for the electron density as a 
function of the fiare peak temperature alone, i.e., UeiTp) oc T^^ (Eq. 34), as it can be seen 
in Fig. 8 by comparing with the RTV predictions for individual flare loop lengths. This 
scahng law predicts a electron densities of rig ^ 10^ — 10^° cm~^ for nanoflares at Tp = 1 
MK, densities of Ue ~ 10^° — 10^^ cm~^ for typical solar flares at Tp — 10 MK, and densities 
of Ue ^ 10^^ - 10^2 cni-3 for large steUar flares at Tp = 100 MK. 

Is there evidence for such higher electron densities in stellar flare loops? Recent reviews 
(e.g., §10 in Giidel 2004; Ness et al. 2004) of electron density measurements in stellar flare 
loops based on density-sensitive iron fine pairs quote a typical range of ~ 2 x lO^-*^ — 2 x 10^^ 
cm~^. A density of rie ^ 3 x 10^^ cm~^ and a spatial scale of 0.1 stellar radius {L ^ 200 
Mm) was measured in a spatially resolved flare on the eclipsing binary Algol B (Schmitt 
et al. 2003). This range of observed densities (n°^* ^ 2 x 10^^ — 2 x 10^"^ cm~^) brackets 
our theoretically predicted range. We have also to keep in mind that previously measured 
densities from stellar spectroscopy (from He-like triplets of OV II, Ne IX, Mg XI, Si XIII, 
and density-sensitive line ratios) have severe sensitivity limitations for densities above > 10^^ 
cm~^. So our densities predicted by the RTV law are in the same ballpark as the observed 
stellar flare densities. The RTV law seems to be a good prediction tool, and the effects of 
short heating times, which can be a factor of 2 lower than the RTV predicted densities 
(Appendix A; Fig. 9), largely cancel out the effects of short heating scale heigths, which can 
reach up to a factor of ~ 2 higher densities than predicted by the RTV law. 
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5.3. An Universal Scaling Law 

A correlation between the volume emission measure EMp = ri^V and the flare peak 
temperature Tp was extended from solar flares to stellar flares (e.g., Feldman et al. 1995b; 
Stern 1992; Shibata & Yokoyama 1999, 2002). A theoretical attempt was made to explain 
the solar/stellar EMp — Tp correlation with a universal flare model in terms of magnetic 
reconnection by Shibata & Yokoyama (1999, 2002). Using the result of numerical MHD 
simulations of flares conducted in Yokoyama & Shibata (1998), where the flare peak temper- 
ature scales as Tp oc B^/'^Uq ^^"^ l?!"^ (with B the magnetic held strength and no the electron 
density outside the reconnection region), setting the thermal pressure equal to the magnetic 
pressure, lUpkBTp ~ i?^/87r, and assuming Euclidian volume scaling, EMp oc n^L^, they 
arrived at an "universal scaling law" of (their Eq. 5) 

--^-»15Ig)"'(io^)"(io^)"i--'1- 

SO the emission measure scales with a power of EMp oc If we introduce the fractal scahng 
of the flare volume, EMp oc ripL^^ , the universal scahng law of Shibata and Yokoyama (1999) 
takes the following form, 

EMp oc s{4-3D^)^J>W2yW2)0^-2 ^ (53) 

which yields the following coefficients for fractal scahng {Dy ~ 2.0 and Dy — 2.4) and 
Euchdian scaling {Dy — 3), 

2nJ °T5 for Dy = 2.0 




EMp oc= { B-^-'^nl^T^-* for Dy^2A . (53) 
^nl^T^-^ for Dy = 3.0 

Thus, if the magnetic field B and electron density tiq outside the reconnection region do 
not have a systematic scaling with the flare peak temperature Tp, the universal scaling law 
of Shibata and Yokoyama (1999) predicts a scaling law of EMp oc T^^ for the observed 
volume fractal scaling of Dy — 2.4, which is somewhat steeper than our measured values of 
EMp oc rp^-5±o.4 for gtellar flares, or predicted by the RTV law, i.e., EMp oc Testing the 
validity of the universal scaling law of Shibata and Yokoyama (1999) requires also statistics 
on magnetic held strengths B and ambient electron densities no in flares. 



6. CONCLUSIONS 



We compiled directly observed parameters from solar and stellar flares, such as the 
volume peak emission measure EMp, flare peak temperature Tp, flare duration t/, and flare 
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length scale L (the latter only for solar flares). A prominent statistical correlation is found 

between the volume emission measure EMp and flare peak temperature Tp, which scales as 
EMp ^ Tp'^ for both solar and stellar flares. Another recent study demonstrated that the 
flare volume has a fractal scaling, V{L) oc L^-^, rather than the generally used Euclidian 
scaling of V{L) oc L^. Applying the RTV scaling law, combined with the fractal volume 
scaling and the statistical L — Tp correlation L{Tp) oc T^-^, leads directly to a theoretically 
predicted scahng law of EMp oc T^ ^, which explains the observed correlations in both solar 
and stellar flares. 

A second result wc find is an unexplained offset by a factor of about 250 between solar 
and stellar ffares at the same temperature, which is likely due to a selection bias for stellar 
flare events with larger volume flUing factors and larger spatial scales. Interestingly, however, 
this selection bias does not affect the overall EM—T relationship and the lower threshold has 
a similar functional dependence of EMmin Tp"^, probably because the detector sensitivities 
are dropping off with a similar function with higher temperatures. 

A third result is that our model of fractal flare volume scaling provides realistic estimates 
of volume flUing factors, and thus of flare densities. We flnd that the electron densities in 
solar flare loops can be predicted based on our fractal scaling in close agreement to the 
predictions of the RTV law. The agreement of the predicted electron densities with both 
methods agrees always better than an order of magnitude (Fig. 7), although the absolute 
magnitude varies by three orders of magnitude between the smallest nanoflares and the 
largest solar flares, i.e., ^ 10^ — 10^^ cm~^. Since the RTV scaling law (combined with 
the observed L (x Tp correlation) predicts about a linear relationship between the electron 
densities and flare peak temperatures, i.e., Up oc Tp, we expect up to an order of magnitude 
higher electron densities in the largest stellar flares due to the higher temperature than in 
solar flares. 

The determination of correct scaling laws allows us also to infer realistic estimates of 
the (conductive and radiative) flare cooling times, which can be tested from the e-folding 
decay time of individual peaks in (solar and stellar) flare light curves. 

The scaling laws allow us also to eliminate temperature biases in the statistics of the 
total thermal energy of flares, i.e., Et oc UpTpV oc EMpTp/up. Since we flnd that the electron 
density (corrected for a fractal flUing factor) scales approximately as Up oc T^ -^, the thermal 
energy scales approximately as Et oc EMp, and thus the observed total emission measure 
EMp can be used as a good proxy for the thermal flare energy Et- Such unbiased frequency 
distributions of flare energies N[Et) permit us then to determine whether there is more 
energy in large or small flares, an important test for nanoflare heating theories. 
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A. DENSITY COMPARISON OF RTV-LAW WITH HYDRODYNAMIC 

SIMULATIONS 

Here we compare the electron densities predicted by the RTV law (Eq. 19; §4.1) with 
maximum flare densities Umax obtained from hydrodynamic simulations of heated flare loops. 
We make use of a parametric study of radiative hydrodynamic modeling conducted in Tsik- 
lauri et al. (2004), where a flare loop with a fixed half length of L = 55 Mm was heated with 
variable heating functions, specified by 5 different heating rates {Eho =0.60, 3.0, 15.0, 30.0, 
and 60.0 erg cm~^ s"^) applied at the flare loop apex for 4 different (Gaussian) durations 
{theat — 41, 82, 164, 329 s). We show 4 evolutionary curves of the apex temperature Tapex{t) 
as a function of the electron density ne{t) for the lowest and highest heating rate, as well as 
for the shortest and longest heating duration in Fig. 9. We see that the 5 different sets of 
heating rates produce approximately maximum flare temperatures of T^ax ~ 14, 22, 35, 43, 
and 53 MK. All parametric values are also tabulated in Table 1 of Tsiklauri et al. (2004). 
For each of the 20 evolutionary curves we mark the maximum values rimax and T^ax of the 
simulated flare loops with a diamond symbol in Fig. 9 and compare it with the RTV law 
(Eq. 21), which is indicated with a dashed line in Fig. 9. We can consider now the density 
ratio n^ax/n'RTV the flare maximum temperature and find that this ratio is a systematic 
function of the heating time scale theat, but is almost independent of the heating rate Ehq and 
maximum temperature T^ax- This ratio amounts to nmax/f^RTV — 1-09 ±0.06 for th^at — 328 

S, nmax/riRTV = 0.81 ± 0.07 for theat = 164 S, nmax/riRTV = 0.57 ± 0.06 for theat = 82 s, 

nmax/nRTV = 0.35 ± 0.06 for theat = 41 s, so the RTV overpredicts the maximum density 
for small heating time scales, but agrees quite well for longer heating time scales and thus 
provides a good proxi to predict the maximum fiare densities. Applying the RTV law to 
scaling laws of fiare parameters, one has to correct for a numerical factor for shorter heating 
times, but this factor is largely independent of the flare temperature. 
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Table 1. Peak emission measure and temperature values in 8 stellar flares. 
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Fig. 1. — Scatterplots of flare peak emission measure EMp{Tp), peak luminosity Lx{Tp), flare 
duration Tf{Tp), and total X-ray radiated energy (Tp) versus the flare peak temperature 
Tp of 68 stellar flares listed in Table 4 of Giidel (2004). Linear regression fits are shown with 
thin lines [correlating y(x) and x(y)] and with thick lines (ordinary least squares bisector 
method), and the linear regression coefficient (RC) is indicated. 
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Fig. 2. — Evolutionary phase diagrams of temperature T(t) versus emission measure 
EM{t) in eight flare events from the foUowing stars: HR1099 [or V711 Tau] (Stern 1996), 
1EQ1839.6+8002 (Pan et aL 1997), II Peg (Doyle et al. 1991), UX Ari (Tsuru et al. 1989), 
Algol (Stern et al. 1992), AB Dor (Maggio et al. 2000), and Proxima Centauris (Reale et 
al. 2003). The straight hne represents the statistical scahng law EM^ oc T^-^ obtained from 
Giidel (2004), which predicts a flare temperature Ts at the maximum emission measure EMp. 
For comparison, also the measured flare peak temperature Tp and the maximum temperature 
Tm are indicated, listed together in Table 1. The ellipses shown in 5 cases represent the 90% 
confidence intervals in two-parameter space, while the temperature error bars in the other 3 
cases are SDecified in the Quoted oaoers. 
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Fig. 3. — Compilation of flare peak emission measure EMp versus flare peak temperature 
Tp measurements in solar and stellar flares. Both the solar and stellar data sets fit a similar 
statistical correlation of EMp oc T^ -^, but the emissison measures of stellar flares are about 
2 orders of magnitude higher at the same temperature. The linear regression fits are shown 
with thick solid fines, while the Icr-ranges that include 67% of the datapoints are indicaded 
with dashed lines. 
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Fig. 4. — Compilation of flare durations rj versus flare peak temperature Tp measurements in 
solar and stellar flares. The combined data sets fit a statistical correlation of tj oc 7"° ■^^^^■'^^ 
(thick solid line), including 67% of the datapoints within a factor of ~ 3 (dashed lines). 
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Fig. 5. — Compilation of flare length scales L versus flare peak temperature Tp measurements 
in solar and stellar flares. The combined data sets fit a statistical correlation of L oc 
(thick solid hne) , including 67% of the datapoints within a factor of f^i 3 (dashed hnes) . 
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Fig. 6. — Conductive cooling times Tc{L,Tp) (top left panel), radiative Tc{L,Tp) cooling 
times (top right panel), combined cooling times l/r — 1/tc + l/r^ for the RTV scahng law 
(bottom left panel) and Serio's scaling law (bottom right panel) are shown as a function of 
the observed flare duration Tf, calculated with the scaling law using the observed flare peak 
temperatures Tp and flare length scales L, and a ratio of L/s// = 3 of the loop half length L 
to the heating scale heigth sh, the maximum limit in Serio's scaling law. 
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Fig. 7. — Comparison of the electron densities determined in nanoflare and flare datasets 
where the flare peak temperature Tp, the total volume emission measure EMp, and the 
flare length scale L = A^l'^ has been observed, using the theoretical RTV scaling law re- 
lationship = Co{Tp/L) (Eq. 50) (x-axis), and the observational relationship n°'^^ = 
^jEMpI {qvLl{L/Lo)^y) (Eq. 49) (y-axis) with Dy = 2.4 and a volume filhng factor range 
of qv — 0.03 — 0.08 at Lq = 25 Mm (vertical error bars). Both methods predict a range of 
He ~ 10^ — 10^^ cm~^ and are consistent within less than an order of magnitude. 
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Fig. 8. — Electron densities predicted by the RTV law as a function of the flare peak 
temperature, for those datasets with measurements of the flare peak temperature Tp and 
length scale L (identical datasets as used in Fig. 5). The straight line is the statistical 
prediction Up = 10^^-^{Tp/10 MK)^-^ (Eq. 34). The range of electron densities (iQH-SiO-e 
cm~^) inferred for stellar flares according to the compilation of Giidel (2004) is indicated 
with a large cross. 
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Fig. 9. — Evolutionary curves of Tapex{t) versus the electron density ne{t) in 20 hydrodynamic 
simulations of ID flare loops from Tsiklauri et al. (2004). The parameters are given in 
Appendix A. All simulations apply to a flare loop with a length of L = 55 Mm. We show 4 
evolutionary curves for the extremal heating rates Ehq = 0.6 and 60.0 erg cm~^ s~^ and for 
the extremal heating time scales t^eat = 41 and 329 s. The locations of T„iax and Umax are 
indicated with diamond symbols. The RTV law is indicated with a dashed line. Note the 
ratios of the maximum flare densities rimax to the values urtv predicted by the RTV law, 
listed at the top of the Figure. 



